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ABSTRACT 



Experimental results using roast turkey meat were 
obtained substantiating the uniformly retreating ice front 
model used in freeze-drying. 

A one-dimensional model for rehydration was developed, 
based on flow in capillary tubes. Experimental results 
using roast turkey meat compared favorably with model pre- 
diction for an effective capillary diameter of 0.0010 cm. 

The results indicate a rapid rehydration following the model 
behavior predicted, followed by a much slower rate of 
rehydration. Maximum rehydration of the samples tested was 
to approximately two thirds of the original water content, 
eighty percent of which occurred during the first, rapid 
rehydration period. 

The experiments also confirmed the model prediction of 
increased rehydration with increasing temperature. 
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I. INTRODUCTION 



Freeze-drying is a process in which the water frozen 
within a substance is removed by sublimation. In order to 
ensure that the water passes directly from a solid phase to 
a vapor phase without passing through a liquid phase the 
temperature of the substance is maintained below the triple 
point temperature of water (32.02°F or 491.69°R). 

There are several advantages to dehydration through 
sublimation rather than through evaporation. In the field 
of food preservation freeze-drying is superior to every 
other drying method in producing the highest quality 
product. In contrast to the shrunken, shrivelled product 
of evaporation, the freeze-dried specimen essentially 
retains its previous size and cellular structure. This is 
because of the structural rigidity provided by the frozen 
state of the material. This porous, non-shrunken structure 
will easily and almost completely re-hydrate. 

The* low processing temperature involved in freeze-drying 
helps to minimize the extent of various degradative chemical 
reactions which often occur during drying processes such as 
non-enzymatic browning protein denaturation and enzymatic 
reaction [Ref. 1]. 

The relative absence of liquid water also helps to 
minimize degradatiys reactions and in addition prevents 
relocation of any soluble species during drying. 
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Another advantage of freeze-drying is that the local 
area undergoes a rapid transition from a completely hydrated 
to a nearly completely dehydrated condition. This also 
contributes to the minimizing of the degradative reactions 
already mentioned. 

Freeze-drying has been used for the preservation of 
biological specimens for many years. Freeze-drying was used 
for the preservation of blood-plasma during World War II 
[Ref. 2]. At the present time the largest applications of 
freeze-drying are in the field of food preservation. Many 
other applications have been suggested and explored, 
including ... "the dehydration of radioactive wastes, the 
preparation of porous catalysts, the stabilization of free 
radicals, and the sublimation of non-aqueous solvents to 
allow operation in different media at different temperatures." 
[Ref. 1]. Noyes [Ref. 4] provides a summary of the progress 
and developments in freeze-drying by many organizations 'and 
corporations in several countries. 

Freeze-drying has been described in general in works by 
Cotson and Smith [Ref. 51 , Rey [Ref. 3] and Van Arsdel 
[Ref. 6], Comprehensive reviews were published by Harper 
and Tappel [Ref. 7] in 1957 and by Burke and Decaureau 
[Ref. 8] in 1964. In 1963 the United States Department of 
Agriculture published a list of 6 33 freeze-drying references 
which had been published prior to 1963 [Ref. 9]. A summary 
of the work on the drying of solids in the Soviet Union was 
published by Fulford [Ref. 10] in 1969 . 
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II. PHYSICAL MECHANISM OF 



FREEZE-DRYING 



An analytical evaluation of the mechanism of freeze- 
drying is beneficial from the viewpoint of determining the 
parameters that affect freeze-drying rates and how changes 
in these pertinent parameters result in changes in the 
rates. These considerations are important because of the 
requirements of both science and industry to provide an 
economical, consistent, high-quality product. 

A. NATURE OF THE FROZEN FRONT 

Most qualitative analyses of the mechanism of freeze- 
drying have assumed that a discrete surface divides the 
frozen hydrated region from the dehydrated region. The 
assumption is that initially a very thin layer at the surface 
nearly completely sublimates and that the frozen front then 
moves inward, leaving a dehydrated solid material. This 
idealization had been universally assumed and accepted. 
Recently, however, discussion has occurred challenging the 
validity of this assumption. 

1. Arguments Against a Sharp Frozen Front 

Meffert [Ref. 11] concluded that a significant 
percentage of the initial water content was not removed by 
the frozen front, but was left behind and finally removed 
by secondary desorption. 

Bralsford [Ref. 12] inferred that if any liquid 
were present a liquid diffusion effect could exist. He 
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concluded that this might undermine the validity of the 
sharp frozen front model. 

Brajnikov et al. [Ref. 13] performed several 
elaborate experiments and concluded that a transition zone 
existed in the sublimation front. 

Luikov [Ref. 14] reports in great detail the 
behavior of a sublimation front. These studies can account 
for a transition zone of, at the most, one millimeter in 
thickness . 

In each of the experiments previously cited there 
were undermining physical factors that precluded the 
definitive conclusion of a diffuse sublimation front. In 
several experiments sample layer thickness and thermocouple 
thickness were of the same order of magnitude as the thick- 
ness of the reported transition zone. In other cases 
melting occurred, because of improper temperature control 
or due to the presence of artifical salts which had been 
introduced for experimental purposes. 

2 . Arguments For a Sharp Frozen Front 

Many studies relate the physical visual observation 
of a sharp discrete frozen front separating the frozen fully 
hydrated portion of the sample from the almost fully 
dehydrated uniformly discolored portion. This was observed 
and reported by Harding [Ref. 15] in beef, Margaritls and 
King [Ref. 16] in turkey meat and Beke [Ref. 17] in pork. 
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In addition to physically observing this phenomenon 
Beke shows photographs which reveal a distinct frozen front 
separating the dry zone from the darker frozen zone. 
Photographs are also provided by Clark [Ref. 18]. These 
also show a sharp demarkation between the frozen and 
dehydrated regions of the specimen. 

Hatcher [Ref. 19] used gamma ray attenuation to 
measure the moisture profile across his sample. Visual 
inspection revealed a planar ice-front. In addition the 
Ice front was observed to be sharper than the 3/16 inch 
diameter of the gamma ray beam. The gamma ray count indi- 
cated no residual moisture once the sublimation front had 
passed. Of particular interest are Hatcher’s reports of 
temperature histories. during the freeze-drying of beef slabs. 
Measurements were made with thermocouples embedded in the 
body of the samples. Figure 1 is a plot of temperatures 
versus time for various thermocouple positions. Figure 2 
is a plot of the same temperatures versus thermocouple 
distance from the sample surface for various times. From 
these plots it is observed that a sharp transition occurs 
from the frozen zone of uniform low temperature to the dry 
zone with a nearly linear temperature gradient. 

Similar results have been observed by Speiss 
[Ref. 20] in the freeze-drying of potato starch and egg 
white and by Sandall et al. [Ref. 21] in the freeze-drying 
of turkey meat. 
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Figure 1. Temperature versus time at different thermocouple locations 
during the freeze-drying of beef. Hatcher [Ref. 19]. 
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Figure 2. Temperature versus position at different times during the 
freeze-drying of beef. Hatcher [Ref. 19 J* 



B. CONCLUSIONS 



King [Ref. 1] concludes that there are factors such as 
structural irregularities and molecular and crystallitic 
phenomenon which may cause a diffuse sublimation front. 

He points out, however, that they operate over short dis- 
tances. It is acknowledged that, as the temperature inside 
the material approaches the melting temperature, the frozen 
front will become more diffuse. However, in the temperature 
ranges which result in a high-quality freeze-dried product 
the supposition of a sharp sublimation front is valid and 
in fact more accurate than many other assumptions made in 
the design and analysis of freeze-drying. 
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III. THEORY OF FREEZE-DRYING 



Figure 3 [Ref.'l] shows the several factors which 
interact to govern the freeze-drying rates. A heat source 
is necessary to provide the heat of sublimation of 1200 BTU 
per pound of water. A moisture sink is required to remove 
the generated water vapor. The driving force for the heat 
transfer is the temperature difference between the heat 
source and the sublimation front. The driving force for 
the mass transfer is the partial pressure difference between 
the partial pressure of the water vapor in equilibrium with 
the sublimation front and the partial pressure of water vapor 
in equilibrium with the moisture sink. 

A. THE FREEZE-DRYING MODEL 

The model most often used in an analytical study of 
freeze-drying is the Uniformly Retreating Ice Front Model 
(URIF) shown in Figure 4 [Ref. 22]. It is assumed that the 
heat ransfer and the mass transfer occur at a pseudo 
steady-state across the dry layer. It is further assumed 
that .the linear temperature profile across the dry surface 
is not distorted by the outflow of water vapor. Sandall 
[Ref. 21] and Sandall and others [Ref. 22] made theoretical 
studies of these assumptions and determined that the maximum 
error introduced by either of these assumptions was four 
percent of the predicted rate of heat or mass transfer. 
Finally it is assumed that complete drying occurs during 
the retreat of the sublimation front. Sandall and others 
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Figure 3. Factors influencing rates of freeze-drying. King [Ref. 1]. 




Figure 4. Model of a uniformly 
retreating ice front 
for the freeze-drying 
process. King [Ref. 1], 
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[Ref. 22] determined that a maximum error of three percent 
resulted from this assumption. 

For these conditions, the expression for heat transfer 
is [Ref. 22] 

q - h(T e - T s > - £ ( Ts . T f ) (1) 



where the heat transfer to the substance by convection is 
equal to the heat transfer through the substance by conduc- 
tion. The convection heat transfer is from the external 
heat source to the surface of the specimen. The conduction 
heat transfer is from the surface of the specimen to the 
sublimation front. 

The surface temperature, T , can be eliminated from 
Equation (1) resulting in 



q 



k 

! + 





( 2 ) 



Similarly the mass flux has the following expression 
[Ref. 22] 



MU = _£ / _ \ _ D]_ p fw ~ p 

Aw R ip sw ew AL 



sw 



( 3 ) 



where the mass transfer external to the specimen is equal 
to the mass transfer internal to the specimen. The mass 
transfer internal to the specimen is from the sublimation 



22 



front to the specimen surface. The mass transfer external 
to the specimen is from the specimen surface to the moisture 
sink. 



Note that T is the average of the absolute temperature 
of the specimen surface and the frozen surface. The 
assumption that T is constant will result in no more than 
two percent error using the extremes of temperatures 
experienced in this experiment. 

The surface partial pressure can be eliminated from 
Equation (3), resulting in 



Since it is assumed that complete drying occurs behind 
the sublimation front, the fraction of water remaining in 
the substance, x, is equal to the fraction of the substance 
still frozen. Thus: 




(4) 



k + AL 
c 




X 



L 

2 



Solving for AL yields 



AL = | (1 - x) 



(5) 
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Substituting Equation (5) into Equation (2) yields 



' ' | + (1-x) £ 



L ^ T e " V 



( 6 ) 



Substituting Equation (5) into Equation (4) yields 



„ w _ D' (p fw “ p ew ) 

N A^ - (1 . x) | 



RT 



(7) 



Assuming that the temperatures in the specimen are not 
changing rapidly it follows that 



9 ■ iH s "w N A 



(8) 



where the heat flux from the specimen is equal to the 
product of the molar flux times the molecular weight times 
the latent heat of sublimation. 

Multiplying Equation (7) by H A and dividing Equation 
(7) by Equation (6) yields 

p fw ” p ew k RT 1 

T - T- D' AH 

e f w s 



F- + I 

c 

- + (1-x) - 

h. k - L X) 2 



(9) 



This is known as the "wet-bulb" relationship for freeze- 
drying [Ref. 22]. 

Evaluating this "wet-bulb" relationship for conditions 
at the surface of the sample leads to [Ref. 21] 
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( 10 ) 



p fw “ p sw k RT 1 
T s " T f ' D ' *w AH s 



From Equation (10) it Is concluded that fixing T and 
therefore P sw > results in a decreasing p^ w with an 
increasing T^. However, at saturation the partial pressure 
increases with increasing temperatures. Thus, the tempera- 
ture and partial pressure must remain constant to satisfy 
Equation (10) and the natural phenomenon involved. 

Applying a similar argument to Equation (9) it is 

concluded that a constant T and p „ will also exist for a 

s sw 

constant T g and p gw if h and k c are sufficiently large when 
compared to the internal coefficients. 

The drying rate can then be related to the heat and 
mass transfer rate equations [Ref. 22]. 



Heat : 



L ^s , dX\ 

2 V w ^ d0 ' ~ 



F + I 



(T e - T f ) 



( 11 ) 



Mass 



L _1_ / dXx 

2 V w de 



D' . ,, s L 
^ + (1 “ X) 2 



S (p f“ ' p « ) (12) 



These equations, (11) and (12), can now be rearranged 
into forms more useful for testing the validity of this 
model and for attaining the transport coefficients from 
drying rate data [Ref. 22] 



25 
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(13) ' 



and 



(1-x) 



' vM 
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(14) 



Thus it is anticipated that a plot of (1-x) versus 



f 

< 

c 





} 



would be linear with the slope dependent 



on the thermal conductivity and the intercept dependent on 
the external heat transfer. 

, ^ f Pf w “ p 

Similarly, a plot of (1-x) versus < — ■ ■ - - 

l vTdf) 

should be linear with the slope dependent on the internal 
mass transfer coefficient and the intercept dependent on 
the external mass transfer coefficient. Because of the two 
variables in the second term of Equations (13) and (l4) it 
is necessary to solve for the slope first, solve for k and 
D’ respectively, then evaluate the intercept and solve for 
h and k with k and D' already known. 

V 

In addition to the assumptions and simplifications noted 
in Section III. A. it is assumed that internal energy absorbed 
by the dry layer to increase its 'temperature as freeze-drying 
progresses is negligible. Sandall [Ref. 21] evaluated this 
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assumption and found that the maximum error for the most 
extreme conditions of freeze-drying was four percent. 

Similarly it is concluded that the energy required to 
heat the vapor as it flows through the dry layer is 
negligible . 

B. TRANSPORT PROPERTIES 

1. Thermal Conductivity 

Figure 5 [Ref. 21] is a plot of the thermal conduc- 
tivity of freeze-dried turkey meat versus pressure. These 
measurements were made in a steady heat flux, thermopile 
apparatus. At very low pressure where the mean free path 
of the vapor is much greater than the pore diameter of the 
specimen it is observed that the vapor provides no signi- 
ficant contribution to the overall thermal conductivity. 

As pressure increases the mean free path of the vapor 
decreases until it is on the order of magnitude of the pore 
diameter. Continued increase in pressure results in a 
linear decrease in the mean free path and a linear increase 
in thermal conductivity. At high pressures the mean free 
path of the vapor is very small compared to pore diameter 
and the effective thermal conductivity does not vary with 
pressure. 

2. Mass Transfer 

Figure 6 [Ref. 21] is a plot of the internal mass 
transfer coefficient of turkey meat versus total pressure. 
The average Knudsen and Bulk Diffusion coefficients of 
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Figure 5. Thermal conductivities of freeze-dried turkey breast in dry air 
Sandall and others [Ref. 22] 
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Figure 6. Internal mass transfer coefficient for freeze dried 
turkey breast in N ? . Sandall [Ref. 21]. 



turkey meat were measured by Gunn [Ref. 22] and plotted. It 
is observed that at the lower pressures the flux is maximum. 
As pressure decreases, the mean free path of the vapor 
increases with the result that the flux varies inversely 
with pressure. 

C. THE LIMITS 

It is necessary to acknowledge that limits exist for both 
the heat transfer driving force and the mass transfer driving 
force. There is an upper limit to the outer surface 
temperature, T , beyond which, thermal degradation of the 
dried product will occur. The upper limit for most meats is 
about 60°C. [Ref. 6]. The second limit is the upper limit 
to the frozen zone temperature, T~, above which the undried 
portion would no longer be frozen. Of interest, this is 
not the melting point of water but normally several degrees 
lower. Sandall, [Ref. 21] for example, found that for 
turkey meat the maximum ice temperature was -5°C, primarily 
because of dissolved solids. 

It is observed in Figures 5 and 6 that at very low 
pressures k is small and D’ is large. Applying these 
factors to Equation (10) and increasing T to increase the 
drying rate, it is observed that the outer surface tempera- 
ture, T , will reach its limit before the ice temperature, 

T^.. Thus the drying rate is more effectively increased by 
increasing k than by increasing D'. Thus freeze-drying at 
low pressure is heat-transfer controlled. 
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Again from Figures 5 and 6 it is observed that at 
relatively high pressures, k is large and D* is small. 
Applying to Equation (10) and increasing T to increase 
the drying rate, it is observed that the ice temperature, 

T^., will reach its limit before the outer surface temperature 
T . Thus the drying rate is more effectively increased by 

o 

increasing D' than by increasing k. Thus high pressure 
freeze-drying is mass-transfer controlled. 
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IV. THEORY OF REHYDRATION 



Relatively little study or research has been undertaken 
in the area of rehydration of freeze-dried foods. 

King, Lam and Sandall [Ref. 24] determined that 
conditions which resulted in a lack of shrinkage during 
the freeze-drying of turkey also resulted in improved 
rehydration. These findings were confirmed by Beke 
[Ref. 17] as also applicable to pork. 

Karel [Ref. 25] observed a decrease in the rehydration 
of haddock fillets freeze-dried with an increase in the 
frozen zone and outer surface temperatures. Karel also 
observed that for freeze-dried salmon the obverse was true; 
improved rehydration occurred for samples processed with 
higher frozen zone and outer surface temperatures. 

Goldblith [Ref. 3] determined that for freeze-dried 
shrimp higher frozen zones and outer surface temperatures 
resulted in reduced rehydration. 

Karel [Ref. 25], Goldblith [Ref. 3] and Cotson and 
Smith [Ref. 5] observed that raw meat and fish were much 
superior in quality when rehydrated with cold water than 
with hot water over a range from 40°F to 212°F. The samples 
rehydrated with hot water had a tough fibrous texture and 
a "twice-cooked" flavor. This result was thought to be 
because of the protein denaturation resulting from the 
higher rehydration temperatures. 
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Hamm [Ref. 2 6] also observed the decrease in the water 
holding capacity of freeze-dried beef with increasing 
rehydration water temperature and observed that a significant 
decrease in rehydration began at 86°C. 

A. THE REHYDRATION MODEL 

Margaritis and King [Ref. 16] developed a structural 
model of a single fiber based on their microscopic measure- 
ments and supplemented by histological information given by 
Maximow and Bloom [Ref. 27]- This model is illustrated in 
Figure 7. 

The wet fiber is contained by the cell wall, the 
sarcolemma, a membrane of one micron thickness. The fiber 
is made up of myofibrils one micron in diameter with a 
one-tenth micron separation distance. The fiber itself has 
a diameter of 36 microns. In the unfrozen state the 
myofibrils are distributed evenly throughout the fiber. 

Figure 7(a). 

During freezing the water moves to the center of the 
fiber, forming an ice core and forcing the myofibrils 
toward the sarcolemma. 

When frozen the fiber has a single ice column of 24 
microns diameter at the center of the column. The myofibrils 
are shrunken and closely packed into the 5 micron annulus 
lying between the sarcolemma and ice core. Figure 7(b). 

The freeze-dried fiber has the same dimensions as the 



33 





Sorcolemmq Wall (I/O 
Myofibrils 
Ice Core 
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(c) Freeze - Dried Fiber 



Figure 7. Cross-section of a model fiber. 

Margaritis and King [Ref. 16]. 
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frozen fiber with the ice column removed from the center 
of the fiber. Figure 7(c). 

Koonz and Ramsbottom [Ref. 28] studied the formation of 
crystals in poultry tissue using a wide variety of freezing 

l 

rates. Their conclusions verify the validity of this model 
of a single ice column in each fiber. 

B. CAPILLARITY AND SURFACE TENSION 

The above rehydration model consisting of a single ice 
column in each fiber suggests a model involving capillarity 
and surface tension as illustrated in Figure 8. 

Moore [Ref. 29] and Barrow [Ref. 30] developed the 
equation for surface tension. This involves a static 
balance between the force of surface tension in the upward 
direction and the force of gravity in the downward 
direction. 

irdy cos <f> = ( ) p S h , ' ( 18 ) 

Solving for surface tension yields 



Y 



d _Pghj_ 
4 cos <j> 



(19) 



For water the contact angle, <J>, is zero [Ref. 29] 
resulting in 



Y = 



d . 
If P S h 



( 20 ) 
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FIGURE 8. Capillary Model 
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Equation (20) is an expression for surface tension under 



steady state conditions which will describe the ultimate 
degree of rehydration. 

It is also desirable to formulate an expression des- 
cribing the rate of rehydration. In this case the difference 
of the terms for force of surface tension and force of 
gravity is taken in order to obtain the net upward force. 



Streeter [Ref. 31] gives the following expression for 
pressure drop for laminar flow in a pipe 

AP s 32 U AVG y ( 

h * g a 2 
& c 

The net upward force due to the differential pressure 
is that differential pressure multiplied by the area 



Combining Equations (22) and (23) yields the following 
expression for F z 





( 21 ) 




(23) 



F 




(24) 



z 
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And combining Equations (21) and (24) yields 



8ir U AVG y h »' 
So 



ird 4 



= ir dy cos $ - ( — q— )pgh, (25) 



Solving Equation (25) for average velocity yields 



u 



AVG 



dh, = 
dt 



[Trdy cos <p - ( 
Birhyj 
So 



nd 

T- 



) pgh,] 



( 26 ) 



An order of magnitude comparison of the two terms in 
the numerator of Equation (26) indicates that the magnitude 
of the second term is always less than two percent of the 
first term for the samples used in this experiment. 

Thus Equation (26) now becomes 



U AVG 



dh, . d ^c 
dt 8h.y 



(27) 



Rearranging and integrating yields 



/ hj dh, 

0 



h 



2 

i 



dYs c 

~w~ 

dyg c t 



t 

/ dt 

0 



(28) 

(29) 



Recall that Equation (29) resulted from the development 
of a vertical capillary model in which the force of gravity 
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acting against the capillary flow was negligible. This 
means that Equation (29) is valid for a horizontal capillary 
model if the force of gravity acting in the direction of 
capillary flow is also negligible. Numerical evaluation 
further indicates that this assumption is valid for speci- 
mens used in this experiment if rehydration immersion 
depths are not greater than three times the specimen dimen- 
sion in the direction of the capillary. 

Thus this dynamic analysis indicates that a second 
order response is expected. Specifically rehydration 
varies directly with the square root of capillary diameter, 
surface tension and time and inversely with the square 
root of viscosity. 

Furthermore, this analysis indicates that for a given 
group of similar samples of constant diameter, rehydrated 
for the same period of time, the rehydration would vary 
with the square root of the ratio of surface tension over 
viscosity. Since this ratio increases with temperature, 
(Appendix B) it is qualitatively predicted that on the basis 
of temperature change alone, that rehydration increases 
with an increase in temperature. 
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V. EXPERIMENTAL RESULTS 



Experimental rsults were undertaken with several 
objectives. It was desired to experimentally verify or 
refute: the concept of a discrete sublimation front; the 

freeze-drying model and the rehydration model. Furthermore 
it was desired to experimentally evaluate the effect of the 
varying of several parameters on the resulting rehydration. 

A. GENERAL DESCRIPTION OF SYSTEMS 
1. Vacuum System 

Figure 9 illustrates the Vacuum System. 

The vacuum is maintained by a National Research 
Corporation 100 S Rotary Gas Ballast Pump. The five 
horse-power pump has a capacity of 100 cubic feet per 
minute and a blank off pressure of ten to fifteen microns. 

The main vacuum chamber is a modified National Research 
Corporation Model 2555 furnace tank. It is a 40 inch long 
by 40 inch diameter stainless steel cylinder. Access is 
through a door at one end of the cylinder. The door is 
moved by an air-controlled hydraulically-operated hoist. 

The door contains a viewing port. At the opposite end of 
the chamber is a variac controlled light. In the event 
of electrical failure, the main vacuum chamber is isolated 
by the automatic closing of the electrically-controlled 
air-operated high vacuum valve. 
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Figure 9- Vacuum System 




The rotary gas ballast pump removes from the 
main vacuum chamber all non-condensible gases that enter 
the chamber when it is open to the atmosphere. In addition, 
the pump removes condensible gases that enter the chamber 
and also those that sublimate from the sample specimens. 

The pump discharges the non-condensible gases to the 

i 

atmosphere. However the condensible gases, which 
are taken into the pump as vapors, condense during the 
compression portion of the pump cycle. These condensed 
gases contaminate the seal oil and reduce the pump effi- 
ciency and the quality of the vacuum. These condensible 
vapors are removed by the process of gas-ballasting. 
Gas-ballasting admits filtered air to the pump when the 
vacuum chamber is isolated by the rotating vane but before 
compression begins. Air and vapor are then compressed 
together to reach or exceed atmospheric pressure before 
the components of vapor individually reach the pressure at 
which they would condense. 

2. Refrigeration System 

Figure 10 illustrates the refrigeration system. 

The main chamber is cooled by a one-third horsepower 
Tecumseh Products Company refrigeration plant charged with 
Freon 502. The refrigerant cools the main vacuum chamber by 
passing through stainless steel coils located in a 36” 
diameter cylindrical shell inside the chamber. 

A condenser is attached to the main vacuum chamber. 
The condenser is cooled by a one-quarter horsepower Bendix- 
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Figure 10. Refrigeration System 
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Westinghouse Corporation refrigeration plant charged with 
Freon 502. The refrigerant cools the condenser by passing 
through three sets of concentric copper coils. 

3. Electrical System 

All power provided to the freeze-drying plant is 
110 volt, 60 cycle A.C. 

In addition to the electrical components already 
mentioned there is an electrically-controlled air-operated 
fore line' valve . This valve closes automatically on loss 
of electrical power. 

4. Electronic System 

Figure 11 illustrates the electronics system. 

Pressure indication is provided by a Veeco Instru- 
ments, Incorporated, Thermocouple Gauge Control. 

Temperature indication is provided by a Leeds and 
Northrup millivolt potentiometer. 

Figure 12 illustrates the electronics associated 
with the weight change indication device (WCID). The 
WCID consists of the specimen tray suspended by a linear 
spring. The bottom of the tray is attached to the movable 
core of a linear core transducer. The transducer is 
excited by an audio oscillator. The transducer output is 
read on a Hewlett-Packard 3^00 A voltmeter and a Mosely 
Autograf Model 7100 B strip chart recorder. 

For a detailed description of the weight-change- 
indication device see Appendix C. 




Figure 11. Electronic System 
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Figure 12. Weight Change Indication Device 
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5 . System Operation 



Prior to an actual freeze-drying evolution, all 
system components are operated until temperatures and 
pressures are at equilibrium. 

The system is capable of maintaining a vacuum of 
250 microns, a chamber temperature of -15°C and a condenser 
temperature of -30°C. 

Next the frozen specimen which is at approximately 
-25°C is placed in the chamber and freeze-drying commences. 
The temperature difference driving force is the temperature 
difference between the heat source and the surface of the 
specimen. The water vapor partial pressure difference 
driving force is the difference between the partial pressure 
of water vapor in equilibrium with the sublimation front 
and the partial pressure of water vapor in equilibrium with 
the moisture sink. In this experiment the heat source is 
the inside chamber walls which are heated by the atmosphere. 
The moisture sink is the condenser. As discussed previously, 
the vacuum pump is capable of removing water vapor from the 
system. However, the primary purpose of the vacuum pump is 
to decrease the chamber pressure. This causes an increase 
in the mean free path of the vapor. As discussed in 
Section III.B.2 this results in an increased mass flow rate. 

Figure 13 illustrates the inter-relationship of the 
various systems. 
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Figure 13. Freeze-Drying Plant 




B. EXPERIMENTAL PROCEDURE 



1. Preparation of Sample 

Turkey breast meat was the specimen selected for 
this experiment. 

The turkey was cooked for 110 minutes at 400°F. 

The turkey was then cut into sample size pieces. 

Two sample sizes were used. Turkey roast slices, approxi- 
mately one centimeter thick and ten centimeters in diameter 
were used for the one-dimensional studies (discrete 
sublimation-front analysis and evaluation of the models for 
freeze drying and rehydration). 

Turkey roast cubes approximately two centimeters on 
a side were used for other rehydration studies. 

The samples were then frozen in a freezer in which 
temperature was maintained between -25°F and -15°F. 

2. Processing of Sample 

The samples were then freeze-dried. During freeze- 
drying, measurements were taken of specimen temperatures, 
specimen weight change and time. 

After completion of freeze-drying, rehydration 
experiments were conducted on the cube samples. It was 
desired to determine the effect of varying several parameters 
on rehydration. 

One of the variables was the extent of freeze-drying. 
Some samples were freeze-dried only until there was no 
further weight-loss. Other samples were left in the 
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freeze-drying chamber for a period of time after the indi- 
cation that weight loss had stopped. 

Another variable Was the time elapsed between the 
completion of freeze-drying and the rehydration process. 

Some samples were rehydrated immediately following the 
completion of freeze-drying. Other samples were exposed to 
the atmosphere for a period of time between freeze-drying 
and rehydration. 

Another variable was the temperature at which 
rehydration occurred. Rehydration temperatures ranged from 
i|0°F to 212°F. 

Another variable was rehydration time. The specimens 
were rehydrated for various periods of time varying from one 
second to three minutes. Following each rehydration the 
specimen was weighed. 

C. EXPERIMENTAL RESULTS 

1. Discrete Frozen Front 

An experiment was conducted in an attempt to verify 
or refute the analytical conclusion of a discrete sublima- 
tion front. 

A turkey slice described in Section V.B. was 
imbedded with four thermocouples placed at even intervals 
in the approximately one centimeter thick slice [Fig. 14]. 

The results are shown in Figure 15 and Figure 1 6 which are 
observed to be similar to Figure 1 and Figure 2. 
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Figure 15. Temperature versus time at different thermocouple locations during the 

. , ,, . , freeze-drying of turkey . . ....... ,,. s . ,, ... », . *. 
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Figure 16. Temperatures versus position at different times during the freeze-drying of turkey 



Figures 15 and 16 are not as dramatic as Figures 
1 and 2 because of several significant differences in the 
experimental procedures. In the experiments conducted by 
Hatcher [Ref. 19], the thermocouple separation distances 
are several times greater than in this experiment. Further- 
more, Hatcher used a variable heat source to maintain a 
constant high surface temperature. In this experiment, the 
constant heat source of the chamber walls resulted in a 
slowly rising surface temperature and thus slowly rising 
ice temperature and dry layer temperature. This slowly 
rising temperature makes it difficult to determine whether 
the observed temperature increase is indeed due to the 
passing of the sublimation front or whether it is simply 
the sample heating up due to an imbalance in the rates of 
energy gain and loss. Adding to the difficulty of discerning 
the passing of the sublimation front are the facts that the 
temperature changes are very small and the temperature 
measuring system is relatively insensitive to small tempera- 
ture changes which in this case are very significant. 

Another problem in the area of temperatures is the 
disparity between the indicated temperature and the actual 
temperature. Note that the temperature profile is never 
constant. Instead a gradient exists with the warmer tempera- 
tures toward the surface. This is to be expected, but only 
after the sublimation front has passed. It is felt that the 
thermocouple lead is conducting heat to the thermocouple 
bead causing it to indicate a higher temperature than 
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actually exists. Furthermore, any opening alongside the 
thermocouple lead would permit the entrance of vapor and 
warm air or would be a surface through which freeze-drying 
could occur resulting in a higher temperature at the 
thermocouple bead. 

2 . The Freeze-Drying Model — The Uniformly 
Retreating Ice Front (URIF) 

An experiment was conducted in order to verify 'or 

refute the analytical freeze-drying model. 

A turkey slice described in Section V.B. was placed 

on the weight-change indication device and freeze-dried. 

It is anticipated from Equation (13) that a plot 
T e ' T f 

of ( 1— x ) versus ^ would be linear. Since T and T~ 

( \ ^ ^ 

’ d6 

are essentially constant during the removal of the first 
85# of the water. Equation (13) may be integrated with a 
constant (T & - T^) resulting in the following form: 



(1-x) 



8kV w 

AH L 2 
s 




(30) 



0 

Thus a plot of (1-x) versus = should result in a 

straight line for the removal of the first 85 # of the water. 
Such a linear plot would suggest adherence of the experimental 
results to the analytical model. 

0 

Figure 17 is a plot of 1-x versus ^ for the turkey 

-L—X 

slice previously mentioned. The curve is essentially linear 
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for the first 85% of weight loss verifying the validity of 
the analytical model employing the Uniformly Retreating 
Ice Front (URIF) . 

3 . Re hydration 

For the experimental evaluation of the rehydration 
model a Rehydration Number and Rehydration Ratio are defined. 

Let 

M q = mass before freeze-drying 
and 

= mass after freeze-drying 
* 

Then 

= M q - Md = mass of water removed during freeze-drying 

X 

Assume that 

M w = maximum amount of water to be absorbed in rehydration 

Let 

= mass of water absorbed at any time after rehydration 
has begun 
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Figure 17. (1-x) versus 6/(1 -x) during the freeze-drying of turkey 



And 



M = + M = mass of specimen at any time after 

rehydration has begun 

Now define the Rehydration Number, R^, 



R H M 

o 



m , + K 

d w 

M , + M 
d w. 



( 30 ) 



Next define the fraction of rehydration, r, such 



that 



M 



r = 



w 



M. 



( 31 ) 



w 



Thus 



M + r M 
_ d r 1 w x 

H M , + M 

d wx 



( 32 ) 



which reduces to 



1 + r 



^ - 






( 33 ) 



M . 



Define the minimum value of the Rehydration Number, 
that is, when M = M^, as R^, thus 



M. 

TD = d 

n HD M 

o 



( 3 ^) 



58 



Now 




which reduces to 

*H ' R HD + rtl - W (36) 

where (1 - R^) is the fraction of water in the original 
mass before freeze-drying and r(l - R^) is the fraction of 
water in the mass at any time during rehydration. 

Thus, Equation (36) may be written verbally as 



FRACTION OF THE 
ORIGINAL (BEFORE 
FREEZE-DRYING) 
MASS AT ANY TIME 



FRACTION OF 
THE ORIGINAL 
MASS THAT IS 
DRY 



FRACTION (r) 

OF THE FRACTION 
OF THIS ORIGINAL 
MASS THAT WAS 
WATER 



(WATER + DRY 






(WATER - 


SUBSTANCE - 




(DRY SUBSTANCE - 


FUNCTION OF 


FUNCTION OF 
TIME) 




CONSTANT) 


' TIME) 



Referring to Figure 8 and assuming that the capillary 
diameter remains constant. Equation ( 31 ) becomes 



r 



M. 



w 



M. 



w 




( 37 ) 
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Dividing through both sides of Equation (29) by 
L 2 yields 



h 2 = t-VZc* 
L? 4yL 2 



( 38 ) 



Note that for a given geometry, the fraction of rehydration, 
r, is directly proportional to the square root of the ratio 
of (j£). This term, which incorporates all of the "fluid 
properties influence on r, is a function of temperature, 
and is tabulated at the end of Appendix B. Prom this data, 
the theoretical model predicts a significant increase in 
rehydration with increasing temperature. 



Now define the Rehydration Ratio, 




as 



R HD 



M 




M 



M, 

M 

o M 

o 




(39) 



Tables 1-4 and 6-9 contain the experimental and 
reduced data resulting from four different batch-runs. 

In each case the samples were prepared and processed 
as outlined in Section B.l. 

After normal processing Batch #1 (Tray #7) was 
freeze-dried for an additional seven days beyond indication 
that weight removal had stopped. Batch § 1 was then removed 
from the freeze-drier and rehydrated. 
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After normal processing Batch #2 (Tray #8) was 
removed from the freeze-drying process, wrapped In two 
plastic bags and left out ..in the atmosphere. After one 
week of exposure to the atmosphere Batch #2 was rehydrated. 

After normal processing Batch #3 (Yellow Tray) was 
removed from the freeze-dryer and rehydrated immediately. 

Batches #1, #2, and #3 were rehydrated at temperatures 
of 1J0°F, 80°F, 100°F, 120°F , l60°F, and 212°F. The specific 
temperatures used depended upon batch size and temperature 
ranges examined. Because of the large piece-to-piece 
variation ‘in sample specimens it was felt that a minimum of 
three specimen pieces was essential for each set of conditions 
examined. 

After reducing the data for Batches #1, #2, and #3, 
and determining that the most effective rehydration tempera- 
tures were the higher temperatures. Batch '#4 (Green Tray) was 
processed and freeze-dried normally and then rehydrated at 
the higher temperatures; l 60 °F, l80°F, 212°F (simmer) and 
212°F (vigorous boil). 

In all cases rehydration was performed by immersing 
the specimen in water for one, two, three, four, five, ten, 
fifteen, thirty, sixty, one hundred twenty, and one hundred 
eighty seconds. The specimen was weighed after each 
rehydration immersion. 

The significant difference of rehydration ratios 
between batches emphasizes the necessity for proper 
processing in order to attain maximum rehydration. 
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Batch # 3 which was freeze-dried until there was no 



further weight loss and then immediately rehydrated shows 
the highest rehydration ratio. In comparison. Batch #1 , 
which was freeze-dried beyond the time of no further weight 
loss and rehydrated immediately shows a significantly 
reduced rehydration ratio. Thus "over-drying" must be 
avoided if maximum rehydration ratio is to be attained. 

Batch #2 was freeze-dried until there was no further 
weight-loss, but left exposed to the atmosphere prior to 
rehydration. The resulting reduction in rehydration ratio 
is significant. Thus "case hardening" must be avoided if 
maximum rehydration is to be attained. 

Even Batch which was freeze-dried until there 
was no further weight loss and almost immediately rehydrated 
shows a reduced rehydration ratio. Despite extensive 
packaging, indications are that case hardening occurred 
prior to rehydration. 

A comparison of rehydration ratios within batches 
points to some significant conclusions. In the cases of 
improperly processed samples, results from Batches § 1 and 
# 2 indicate that the higher rehydration temperatures result 
in the higher rehydration ratios. Thus the "over-drying" 
and "case-hardening" problems are best relieved by high 
temperature rehydration or "re-cooking". 

A comparison of rehydration ratios within Batch #3 
reveals no conclusive trends. If a conclusion were drawn 
it might be that the middle temperatures afforded the highest 
rehydration ratio. 
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Batch #4 was rehydrated over a range of high tempera- 
tures including boiling; both simmering and vigorous boiling. 
Results do not indicate conclusively an obvious optimum 
rehydration temperature. It does seem indicative that a 
hard boil is preferential to a simmer. The somewhat obvious 
conclusion from Batch #4 that l80°F is the most optimum 
temperature for maximum rehydration ratio is not reinforced 
with the inconclusive results from Batch #3. 

The rehydration model predicted an increased 
rehydration ratio with increased temperature. This predic- 
tion is very strongly verified by the experimental data. 

In the batches of "case-hardened" and "over-cooked" 
specimens the model prediction is the only trend noted. 

For the batches of specimens that were freeze-dried 
only until weight loss stopped and then immediately 
rehydrated, the trend is not as significant. There is an 
increase in rehydration ratio with increasing temperature up 
to the range of 120°F - l60°F. At higher temperatures 
rehydration ratio decreases. This agrees with findings of 
other researchers in this field who found a "protein 
denaturation" effect above l80°F. 

The analytical re hydration model predicts that 
rehydration rate would follow a second order curve for one 
dimensional rehydration. Furthermore, the analysis predicts 
that rehydration occurs in about one second in a sample one 
centimeter long in the capillary direction and that 
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rehydration occurs in between four and five seconds for a 
sample approximately 2.2 cm long in the capillary direction. 

Tables 5 and 10-12 contain the experimental and 
reduced data resulting from four one-dimensional samples. 

Pour turkey meat slices approximately 2.2 centimeters 
thick by ten centimeters in diameter were freeze-dried until 
there was no further weight loss and then immediately 
rehydrated in 120°F water. 

Samples #1 and #3 were rehydrated in a horizontal 
position (capillaries vertical) and Samples § 2 and #4 were 
rehydrated in a vertical position (capillaries horizontal). 

As the model predicted rehydration occurs very 
rapidly. Specifically, the model predicted that for a 
sample of 2.2 centimeters thickness in the capillary 
direction, rehydration should occur in approximately four 
to five seconds. 

Figure 18 is a plot of fraction of rehydration, r, 

versus time. Figure 18a is a plot of rehydration number, s 
M 

, versus time. Figure 18b is a plot of rehydration ratio 
o 

versus time. 

The experimental results suggest that of the total 
amount of water that is eventually rehydrated } a fraction 
(about 80 %) is rehydrated in the first few seconds. An 
obvious break occurs in the experimental curves in Figure 
18 at about r = 0.5, which is approximately 80 % of the r at 
three minutes. Prior to this point, at approximately 
t = 3 seconds, it appears that rehydration follows the model 
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of capillary rehydration. After 3 seconds the rehydration 
rate slows significantly, indicating that the large capillaries 
are filled. The subsequent rehydration would seem to be 
the much slower soaking of the previously dry myofibrils. 

The decreased rehydration rate would seem to be due to the 
smaller effective diameter capillaries. 

Thus for purposes of comparing the experimental 
results with the rehydration model, the experimental results 
should be normalized. This will result in experimental curves 
based on only the water which is rehydrated due to primary 
capillary action, which is the basis for the model. 

Figure 19 is'-the comparison of the theoretical curve 
and the normalized experimental curves. Since the break is 

V 

observed to occur at a r * 0.5, the normalized curves are 
obtained by dividing r by 0.5. It is observed that the 
experimental curves strongly verify the rehydration model. 

In addition they specifically suggest an effective capillary 
diameter of 0.0010 cm. 

From the figures it is observed that the difference 
between rehydration ratios of the specimens rehydrated 
horizontally and those rehydrated vertically is insignificant. 
The percentage difference in rehydration ratios is less than 
four percent. Thus the experimental results verify the 
prediction of the rehydration model that the effect of gravity 
on rehydration was insignificant. 
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Figure 18. Fraction of rehydration, r, as a function of time. (Batch #5) 
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time - seconds 

Figure 18a. Hehydration number, M/M , as a function of time (Batch #5) 
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Figure 18b. Rehydration rate, Rrr/Rp- as a function of time (Batch #5) 
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Figure 19. r/0.5 as a function of time (Batch #5) 



VI. CONCLUSION 



Based on the above studies of both an analytical and 
an experimental nature, several conclusions are made. 

Both analytical and experimental results strongly 
suggest that the nature of the frozen front is a discrete 
and not a diffuse sublimation front, and that the Uniformly 
Retreating, Ice Front (URIF) Model for freeze-drying is 
valid. 

In the areas of rehydration there is less definition 
because of the dearth of work in these areas. Based on the 
analytical and experimental results of this author, 
indications are that the rehydration model developed in 
this paper is a valid approximation. Experimental results 
verified the rehydration model in the prediction of 
rehydration rates, the assumption of the insignificance of 
the effect of gravity and the anticipated change of rehydra- 
tion rate with change in temperature. The experimental 
results compared favorably with the model with an effective 
capillary diameter of 0.0010 cm. 

From the general rehydration studies it is concluded 
that rehydration ratios increase with increasing temperature 
for "over-dried" or "case-hardened" specimens. Temperature 
effects on rehydration ratios are inconclusive for specimens 
that have been freeze-dried without "over-drying" and imme- 
diately rehydrated without "case-hardening." 
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Of significance, it should be noted that the newly 
developed rehydration model accurately predicts rehydration 
rates for the first 80 - 85 % of rehydration. Recall that 
the freeze-drying model also accurately predicted freeze- 
drying rates for the first 80-85% of freeze-drying. 

It is apparent that during the last 15-20% of both 
freeze-drying and rehydration secondary effects become 
significant causing the observed deviations from the model 
predictions . 
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Rehydration Rates by Weight Batch # 1 - Tray #7 
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TABLE 1. Rehydration Rate versus time. Batch #1 
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TABLE 2. Rehydration Weight versus time; Batch 



Rehydration Rates by Weight (GMS) Batch #3 
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TABLE 3. Rehydration Rate versus time; Batch #3 



Rehydration Rates by Weight (GMS) Batch - Green Tray 
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TABLE 4. Rehydration Weight versus time; Batch 



Rehydration Rates by Weight in grams Batch #5 
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TABLE 5. Rehydration Rate versus weight; Batch # 5 
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TABLE 6. Rehy drat Ion Ratio, R T j/R Hn , as a function of time; Batch #1 
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TABLE 7. Rehydration Ratio, R T ,/R Hn , as a function of time; Batch # 2 
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TABLE 8. Rehydration Ratio, R„/R„ n as a function of time. Batch #3 
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TABLE 9. Rehydration Ratio, R H /R Hr ., as a function of time; Batch 
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Fraction of Rehydration (Experimental - « — ) Batch #5 
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Rehydration Number (Experimental - rp ) Batch #5 
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THERMOCOUPLE TEMPERATURE (°C) 



APPENDIX A 
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Above taken from strip chart recorder plot 



APPENDIX B 



SAMPLE CALCULATIONS 

Predicted Rehydratlon (capillary height h, — any direction) 

Assume "effective sarcollemma diameter" is 
one quarter "actual sarcollemma diameter" — 

r 

i.e. d = 0.0006 cm 



.2 d y 

h . = v i s 



h 

L 



4 u L 
g Y d 



For 



d = 0.0006 cm 

2 

y = 70.0 gm/sec 
y = 0.01 gm/cm-sec 
L, = 2.2 cm 



t (sec ) 


h f (cm) 


r 


1 
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0.682 
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0.819 
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Square root of ratio of surface tension over viscosity 
for water. 



Temp Ratio [Ref. 21] Square root of 

°F X of ratio [ 



40 0 

80 0 

100 0 

120 0 

160 0 

180 0 

212 0 



159 


0.399 


273 


0.523 


337 


0.381 


398 


0.632 


527 


0.726 


587 


0.766 


680 
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APPENDIX C 



In order to prevent "over-drying" it is necessary to 
have an end-of-run indication. It is desirable that this 
end-of-run indication be reliable, inexpensive and easily 
recordable. One parameter which can provide end-of-run 
indication is specimen weight, since the evolution of 
freeze-drying is completed when there is no further removal 
of water vapor from the specimen. An analytic balance is 
not satisfactory because it requires continuous local 
operation. A Mettler Scale is prohibitive because of its 
high cost. In fact, a weight-measuring device is not 
necessary. All that is needed is an indication of the 
change in weight. The weight-change-indication-device 
(WCID) was designed to fulfill this requirement. 

Physically the WCID is put together on a chemistry 
stand 30 inches high. At the top of the stand is clamped 
a cross-bar from which hangs a specimen tray. The specimen 
tray is suspended by a spring which allows up and down 
movement of the specimen tray as the specimen weight changes. 
The spring is, of course, linear for the range of the 
specimen weight. At the bottom of the specimen tray is 
attached the movable core of a linear core transducer. The 
static outer cylinder of the transducer is clamped to the 
bottom of the chemistry stand. Excitation to the linear 
core transducer is provided by an audio oscillator. The 
output of the transducer is fed to a voltmeter and a trace 
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recorder. The WCID is extremely versatile due to two of 
its features. By changing the spring, the WCID can be 
adapted to handle extremes of specimen weights. Because of 
this feature the WCID is usable over a greater range of 
weights than either an analytic balance or a Mettler Scale. 
The second feature is the facility for adjusting the gain 
of the amplifier and the null position of the movable core 
transformer. This permits the setting of a wide range of 
sensitivity. 

Figure 20 is an illustration of the WCID. 
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Figure 20. Weight change indicator device 
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